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Abstract 

<r 

. For a periodically perforated structure, for which homogenization takes 

00 ' place in the linear theory of elasticity, the components of the effective 

elasticity tensor depend in general on the geometry of the holes as well 
as on the local elastic properties. These dependencies were shown by 
' Vigdergauz in [3 to be separated in an elementary way for one particular 

class of structures. The original proof of this relation made use of the 
lattice approach to describe periodic functions using complex variables. 
<~i I In this paper we present a proof of the so-called Vigdergauz equations for a 

related class of square symmetric structures. Our proof relies solely on the 
fundamental theorem of real variable calculus. The differences between 
the two mentioned classes of structures are nontrivial which makes our 
result a partial generalization as well. 



> 

0\ • 1 Introduction 

, In [7] Vigdergauz showed an explicit relation between the local and the effective 

' elastic properties for a special class of periodic structures in the linear theory 

Tj" , of elasticity. The result relates the components of the effective elasticity tensor 

' to the local elastic properties as elementary functions involving one geometric 

, constant for each distinct component of the effective tensor. For a subclass of 

the square symmetric planar structures with local elastic moduli {K, G) , the 
relation to the corresponding effective properties {K*, G*, G45) can be written 
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where (^41,^2,^3) are geometric constants that depend solely on the domain 
of a periodicity cell of the structure. 
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The main restrictions defining the class of structures studies in [7] were that 
the periodicity ceh should contain only one hole which is centered strictly inside 
the cell. In addition the boundary of the cell was supposed to be smooth. An 
example of such a structure is shown in Figure [IJa). 

In this paper we will consider the class obtained by relaxing the assumptions 
on the number of holes and the smoothness of the boundary. More precisely, we 
consider the structures with positive definite square symmetric effective tensor, 
which have a cell with Lipschitz continuous boundary. Moreover, we suppose 
that the cell can be transversed by some lines parallel to the coordinate axes. 
The last requirement enables us to give a particularly simple proof of the explicit 
relations, which we call the Vigdergauz equations. 

An example of a structure of the kind we will consider is shown in Figure[ljb). 
This is not of the type studied in [7] both because of the lack of smoothness 
and due to the number of holes in any peridicity cell. Conversely, the structure 
shown in Figure [TJa) is not in the class we will consider because no periodicity 
cell can be transversed by line segments. Hence the classes have nontrivial 
differences and so our result gives a partial generalization of the Vigdergauz 
equations. 

As an illustration we consider now the first two of the above equations. To 
emphasize on the idea of proof we will make formal calculations and leave the 
precise definitions and arguments to the later sections. 

Let C be the domain of one period of a structure in and Y the corre- 
sponding cell, which we assume is a rectangle. We consider the case of linear 
elasticity and use the planar bulk modulus K and shear modulus G as the local 
material parameters. The stress a and the strain e are then linearly related via 
the local version of the Hooke law: 



(1) 



here written in Voigt notation and 712 is the engineering shear strain. The 
matrix on the right hand side is the inverse of the representation of the elasticity 
tensor. We write the two first equations explicitly as 





1/1 1 \ 1/1 1 



1/1 1 \ 1/1 



We assume that the strain components Sy , hence also the stress components 
CTy, are F-periodic, where Y = [OJi] x [0,^2], for some positive real numbers 
li and l2- We also assume global quadratic symmetry. Hence, we have the 
corresponding average relations 

(£11) = \{-^ + ^) + \{-^-^) ' 

= l{-^-^) (-11) + \{-^ + ^) (^22) , 

where (•) represents an average over C or Y. Here K* denotes the effective bulk 
modulus, and G* one of the two effective (transverse) shear moduli. 



2 



We assume that the possible holes in C are such that there exists both 
horizontal and vertical line segments contained in C, joining opposite faces of 
Y . Then we can integrate along a horizontal line in the y-cell, and find the 
value / U22 dx. It is possible to show that this integral is independent of y. 
Similarly, we find, by integrating vertically, that J andy is independent of x. 
Thus, we find that 

(0-22) = J cr22dx, (o-ii) = ^ J cjiidy. 

Due to the periodicity of the deformed structure, J en dx — ui{li,y) — ui(0, y) 
and J £22 dy = U2{x, I2) — U2{x, 0), are independent of x and y. Hence, 

{£22) = j £22 dy, (£11) ^ y £11 da;. (2) 

Let us now assume that the average vertical forces is zero, that is {1J22) = 0. 
Then, integrating along lines without holes we obtain 



{£22) - ^ / = J - ^) ]^ / + i + ^) ^ / dy. 

Let a and h be defined by 

(1 + a)(crii) = ^ J aiidx, = j- J cr22dy. 

Then, by the equations in ^ we obtain that 

(.22) = i(^-^) (-n) + i(^ + ^)M.n). 
Combined with the reduced average relations, 

(-11) = l{-^ + ^) (-11) . (^22) (-11) - 

we find that 

1 1 \ , , I I 

K + g) (^ + '^)=if^+G^' 

K g)^\K^g) ^ K* G*' 
Adding and subtacting these two equations yield, respectively, 

I /I l\ a+b I 



K \K GJ 2 K*' 
I /I l\ a-b 1 



G \K GJ 2 G 
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In summary, we have found the following representations of the Vigdergauz 
constants: 



A2 



Ax 



^ / an + ^ / (722 dy ^ \ 
2(aii) 2 

2(aii) 2 



2' 



The original proof of the above mentioned relations made use of some proper- 
ties of the Weierstrass zeta-function and the Kolosov-Muskhelishvili potentials. 
In this paper we present a proof of the Vigdergauz equations for a class of square 
symmetric structures. Our proof relies solely on the Newton-Leibniz and the 
Green formulas, and a lemma of Michell type. The Michell lemma can be seen 
as a consequence of the theorem of Stokes in the current setting. 

The rest of this paper is organized as follows. In the next section we give 
the precise definition of the class of structures that will be studied and set the 
notation. In Section [3] and |4l we give some representations of the average stress 
and the quasiperiods for the cell problem in the definition of the effective tensor. 
These results are used in the proof of the Vigdergauz equations in Section [5j 



Figure 1: Two periodic structures (a) and (b) with square symmetric effective 
elasticity tensors. Periodicity cells are marked with dashed polygons. 



2 Setting of the problem 

We will consider a periodic structure in the planar linear theory of elasticity. Let 
denote a period of the structure. The associated elasticity tensor is assumed 
to be homogeneous and isotropic. The local constitutive relation between the 
stress a and the strain e will be the standard Hooke law which we write using 
matrix notation as 




(a) 



(b) 




K^G K~G 
K-G K+G 
G 




(3) 
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The two parameters K and G above are the planar bulk modulus and the shear 
modulus, which here are assumed to be positive real numbers so that the matrix 
is positive definite. Some additional assumptions will be made on the domain 

n. 

For a displacement field u = (ui,U2) defined on f2, with values in R^, the 
gradient will be decomposed into its symmetric and antisymmetric parts: 

Vw = e{u) + uj{u). 

Explicitly, the components of these are 

\ f dui duA If dui duj\ 

The stress (j{u) is then defined by equation (jSj. When the displacement field is 
clear from the context, it will be dropped from the strain, the rotation, and the 
stress, writing just s, lu, and a. 

The domain Q is assumed to be a bounded and connected open set in with 
Lipschitz continuous boundary. Moreover, we assume that there correspond 
periods li > such that the periodic extension A of i7, 

A = Intel U ((/ifci,;2fc2)+^^), 
kez^ 

is also connected and open with Lipschitz continuous boundary. The translate 
of the cell (0, li) x (0, ^2) that contains will be denoted by Y. 

Apart from the above standard hypotheses on the periodicity cell, we will 
restrict our study to the following special type of structure. We assume that 
there exist two line segments contained in D, that connects the opposite faces of 
Y. This means that some translates 71 and 72 of (0, ^i) x {a} and {(3} x (0, 12), 
respectively, are subsets of il for some a and (3. Illustrations of a perforated 
cell il inside a cell Y with line segments 71 and 72 are shown in Figure [ijb) 
andl^a). 

We denote by Hp^^.{il) the closed subspace of H^{il) of elements with equal 
traces on the opposite faces of dil f) dY. Let S denote the set of symmetric 
2x2 matrices with real entries. We say that u G H^{^1)'^ is quasiperiodic if 
u — ^a; e Hp^j.{n)^ for some ^ G 5, and ^ is here called a quasiperiod. 

Let / be a function defined on some bounded domain X. The average value 
of / will be denoted hy j^^ f dx = \X\~^ J-^ f dx. Any / on X will be considered 
as a function on R^ by periodic extension followed by extension by zero. The 
relevant example of usage is with g defined on fJ, we will write fyd^^ fo'^ 

The effective elasticity tensor is assumed to be square symmetric and positive 
definite. This means additional restrictions on f2 and which can be explicitly 
stated as follows. The effective tensor, described with parameters K* ^ G*, and 
GI5, defined for ^ e S* by the equation 

/fy fill dx\ (K* + G* K*-G* \ / ai \ 
\jya22dx\ ^\K* ~G* K* + G* U22 , (4) 
\jyO,2dx) \ GtJ\2^,2j 

where a is the unique stress obtained from equation (jS]) with u — ^x € H^^j.{^l)'^ 
being a minimizer of the elastic energy J^e{u) ■ a{u) dx. The assumption of 
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71 




□ □ □□ 



(a) 



(b) 



Figure 2: A cell (a) and its extension together with line segments (b). The 
periodicity cell il in (a) corresponds to the dashed rectangle Y in Figure [TJb). 
The line segments marked with arrows in (a) are examples of 71 and 72. In (b) 
the cell shown in (a) is marked as a colored region and the dashed rectangle 
shows the translated cell Y' . The line segments marked with arrows in (b) are 
the translates 'y[ and 72 of 71 and 72, respectively, from (a) as defined in the 
proof of Lemma [1] 



positive definiteness is equivalent to K*, G*, Gig > 0, which are the scaled 
eigenvalues of the matrix on the right hand side in equation (|4]). The matrix is 
in particular nonsingular. 

The existence of a solution u to the above minimization problem is a direct 
consequence of the Riesz representation theorem and the Korn inequality. The 
displacement field u is by periodicity unique up to translation, which shows that 
the corresponding stress a is unique. Moreover, since the data are smooth, u is 
smooth in some neighborhoods of 71 and 72. A proof of the regularity of u can 
be found in [5] and we refer for contreteness also to Chapter 6]. 

The relation between the local parameter K and G, and the components 
Cijki of the elasticity complianc^il tensor in the standard notation is given by 



Cllll - C2222 -1[k + g)^ - I [k ~ Q J ' '''''' " 4G' 

where the rest of the components are zero, up to the standard symmetrjj^. The 
parameters in relation Q are called the effective planar bulk modulus K* , and 
the effective shear moduli G* and G45 . The relation to the effective compliance 
tensor is, similarly to the above description, given by 

1/1 1\ . 1/1 1\ . 1 



4 \K* G* J ' ''^^ 4 \K* G* ) ' '^'^ 4G|5 

The rest of the components are determined by symmetry as in the local relation. 
For a stress field cr, the divergence and the trace will be denoted as follows: 

/ 9(711 9(712 9(712 , 9(722 \ 

"1^"= "5 ^"5 ' Trcr = (711 + (722. 

\ ox\ 0x2 0x1 0x2 J 



^Thc inverse of the elasticity tensor. 
^For all indices, Cijki = Cjiki = Ciji,^. 
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The outward unit normal of at a Lipschitz continuous boundary will be denoted 
by V, and the normal traction of the stress by av. The stress and the strain are 
sometimes represented by the vectors a = (cth, (T22, cri2) and e = (en, £22, 2ei2), 
as in equation Q and in a sense also in equation (|H). The product is then the 
standard inner product on R^, which coincides with the Frobenius inner product 
of the tensor fields. The relavant example here is the elastic energy above. 

3 Representation of average stress 

In this section we give a representation of the average stress over the periodicity 
cell in terms of the average stress over line segments in the definition of the 
effective elasticity tensor. See FigureUJa) for an illustration of the line segments 
71 and 72. In Figure (^fb), the other line segments used in the proof of the lemma 
below are illustrated. 

Lemma 1. Let u be a quasiperiodic minimizer of J^^e{u)-a{u) dx. Then 



Proof. Since u minimizes the elastic energy, the stress cr = cr(u) satisfies the 
equilibrium equation div a = in the sense of distributions with vanishing 
normal stress ai' at the boundaries of any possible hole in the global structure. 
The quasiperiodicity of u implies the periodicity of e and hence of a. We 
therefore have 



Let j[ be the translate of 71 such that the left endpoint of 7^ is the common 
point of 71 and 72 . Let 72 be the translate of 72 such that its lower endpoint is 
the left endpoint of 7^ . Let Y' be a translate of Y such that 7^ and 72 are the 
lower and left sides of the boundary of Y' n A, respectively. Let VI' =Y' Ci A. 



Since a and diva have components in L^(Sl'), the following Green formula 
holds for all ip e H^{n'f: 



where (•,•) denotes the pairing of H^^^{dQ')^ and its continuous dual. Since 
dist{dY', dA) > and ai^ = on OA it follows from the regularity of u that 
aiy e L^{dn'Y. Hence 




div (7 = in A, 
(Ti^ = on dA. 




G ■ V</3 dx + div a ■ (pdx — {av, (p) , 
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Let 77 S R be such that (p = {xi + 77, 0) vanishes on 72. Then by periodicity 
and the Green formula, 

/ (Jixdx = / (Til da; = / <7v ■ [Ix^Q) dx — li I audx — h / audx. 
Jn Jn' JdY' J-i'^ J72 

A division by |y| yields the first equation in the statement of the lemma. 

The other three equations follow by the same argument when using test 
functions of the types (0,X2 + 77), {x2 + 77,0), and (0,2:1 + rf), respectively, for 
suitable 77. □ 

In the proof of the above lemma used the Green formula which relies on the 
existence of a normal trace operator for elements in L^(r2)" with divergence in 
i^(ri). We refer the reader to [3J Chapter 2] (and 4 ) for its construction. 



4 Representation of quasiperiods 

In this section we give a representation of the quasiperiods of the displacement 
fields in the definition of the effective elasticity tensor in terms the average stress 
over the line segments and the periodicity cell. The idea of the proof is to write 
V = e — a; and then use the Newton-Leibniz formula and that the local elasticity 
tensor is homogeneous by supposition. Hence it is essentially an application of 
the Cesaro formula. 

Let 7( and 72 denote the translations of 71 and 72, respectively, that were 
considered in the proof of Lemma [T] For illustrations see Figure [^Ja) and [2Kb). 

Lemma 2. Let u he a quasiperiodic minimizer of J^e{u) ■ a{u) dx. Then the 
quasiperiod ^ G S* of u satisfies 




Proof. Since u is quasiperiodic we have by the Newton-Leibniz formula and the 
local elastic relation ([3]), 

^iih = ui(b) - ui{a) = [ ^— !-dx= / endx 
Jji 0x1 J^^ 

= ^ a + L + i a - ' 

where a and b denote the endpoints of 71. By Lemma (TJ 

^nh = ^ (/^^ dx + ^^ a22 dx^ + ^ (|^^ .n dx - .2. dx^ . 

A division by h gives the first representation and the second follows by an 
interchange of the indices. 
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Denote by a the common point 71 and 72. With h being the other 
endpoint of ^[ , we have by periodicity, 

^12/1 = M2(&) - "2(0) / TT^dx^ [ ei2dx+ [ UJ2ldx 

= ^ / '^i2dx+ LU2idx = — (T12 dx+ W21 dx. 

An integration by parts gives by periodicity, 

LJ21 dx = liUJ2i{a) - [ ^^^{xi - bi) dx = liUJ2iia) - [ ^^^^xidx. 
r, Jfi 9x1 7^; dxi 

By the local relation ([3]), 

f duJ21 , f (d£i2 d£ii\ 

/ "H — xidx^ / {— — ]xidx 

Jy 0x1 J^; V dxi dX2 I 



1 /■ 5^12 , 1 / 1 1 \ /■ 9(711 , 

x\dx 1" / ~R — ^2 dx 



2G Jy^ dxi 4:\K G J Jy^ dxi 

1 f I 1 \ /■ 9(722 , 

xi dx. 



4 \K GJ Jy dx2 



A division by h yields 



C12 =t^2i(a) + -^/ cri2dx-^l ^a;i 



2G^, 2Gjy dxi 



dx 

1 /• , 1/1 1 \ /• dTrcr , 
= -21(a) + -f^ .,,dx + -^- + -^ -^x, dx, 

where we in the last step used Lemma [I] and that div ct = 0. By periodicity 
—12 (a) = — 12(&) and thus we get by an interchange of the indices, 

62 = -12(a) + -^^ (712 d- + l[jl+G) 4 

By the antisymmetry of uj, taking the arithmetic mean of these two representa- 
tions of ^12 completes the proof. □ 



5 The Vigdergauz equations 

In this section we give a proof of the Vigdergauz equations for the considered 
class of structures. The proof is a straightforward application of the representa- 
tion lemmas of the previous sections to the definition of the effective elasticity 
tensor. The independence of the local elastic properties K and G follow from the 
Michell lemma. We briefly recall the requirements on the periodicity cell. The 

^The point a is also the common point of 7J and 



9 



cell ft is assumed to be such that the effective elasticity tensor is positive definite 
and has square symmetry so that it can be described by three parameters K*, 
G* , and G\^. Moreover, it is assumed that there exist line segments in $7, as 
in Figure llja), connecting the opposite sides without touching the boundary of 
any possible hole in the periodic extension A of O. 

The Michell lemma states that for the pure traction problerrQ in with 
homogeneous and isotropic elasticity tensor on a bounded Lipschitz domain, 
the stress field is independent of the values of the components of the elasticity 
tensor. A reference is given in the end of this section. 

Theorem 1. There exist nonnegative real numbers Ai depending only on Q. 
such that 
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Proof. Let K,G > he given. Let rj E S. Then by the nonsingularity of the 
effective tensor, there exists a quasiperiodic minimizer of J^^e{u) ■a{u) dx with 
unique quasiperiod ^ € S and with average stress, defined by the left hand side of 
equation dU, which is equal to rj represented by (?7ii, ?y22, '712)- The displacement 
field u is unique up to translation. Therefore the stress a is unique. Thus a is 
the unique solution to the corresponding pure traction problem. Hence by the 
Michell lemma, the stress a is independent of both K and G. We proceed in 
three steps. 

Choose a, by specifying r/, such that 

+ (Til da; ^ 0, +- (T22 dx = + (T12 dx = 0. 

Jy Jy Jy 

By equation Q and Lemma [H 

1 , 1 4^11 4 '^ii f 1 l\ 



K* G* j^yUiidx fy'^iidx \K G 
Choose a such that 

+ (Til dx = -j- (T22 dx ^ 0, + (T12 dx = 0. 

Jy Jy Jy 

Then by equation Q and Lemma [H 

1 2^11 1 1 / £1 '^ii d-x 



K* ^11 dx K 2\-^y ail dx 




1 h7 + 7^ • (6) 



^The Neumann problem in which only the normal traction au is specified on the boundary 
and no external body forces are present. In addition, the total force on each hole is required 
to vanish. 
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From the equations ^ and ^ , the two first equations in the statement of the 
lemma follow, since the stress is independent of the local elastic properties. 
Choose a such that 

j- ail '^^ ~ "^22 dx — 0, j- (712 dx ^ 0. 

From equation (U) and Lemma [5] it follows that 



1 1 L^^idx + L^x,dx 1 



G^s G fyaudx \K G 



(7) 



which completes the proof of the separation of variables. 

Left to show is the nonnegativity of the geometric constants. This follows 
directly from the assumption of positive definite effective elasticity tensor. For 
the matrix is symmetric and thus has strictly positive real eigenvalues, so the 
same holds for its inverse. Using the equations ©-([7]) to write the components 
of the effective compliance matrix in terms of Ai, A2, A3, and K, G, we find 
that the eigenvalues are 

_ {K + G)Ai + G _ (K + G)A2 + K _ {K + G) A3 + K 
' " 2KG ' ' " 2KG ' - KG 

for any choice oi K, G > 0. Hence Ai, A2, A^ > 0. □ 

Remark 1. There are natural alternatives to letting a be such that its average 
value over the periodicity cell is some scaled element of the canonical basis of 
R3 as in the proof of Theorem QJ For example, one can choose a such that the 
average stress vector varies over eigenvectors of the effective elasticity matrix in 
equation Q), which are clearly also the eigenvectors of the compliance matrix. 
In this way the deformations are connected to the physical intuition of what 
loads the bulk and shear moduli are supporting and in that sense arrive more 
directly to the Vigdergauz equations as examplified in J^. However, the proof 
given above would not be simplified in any way by making such modification. 

Remark 2. The nonnegativity of the geometric constants Ai, A2, A3 in Theo- 
rem]^ can also be seen as a direct consequence of the Hashin-Shtrikman bounds, 
which for square symmetric tensors can be written 
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where p denotes the positive volume fraction p = \D,\/\Y\ < 1. 

A proof of the Michell lemma for smooth domains can be found in (Sj Chap- 
ter 5]. For Lipschitz domains, the theorem of Stokes can be used. We refer the 
reader to [21 Chapter 2] for an appropriate version of the theorem, which should 
be applied to the simply connected components of a decomposition of il. The 



11 



Hashin-Shtrikman bounds in terms of K*, G*, and that were used in the 
above proof, can be found in W. 

Remark that if some lines at angles 7r/4 and 37r/4 transverse the corre- 
sponding periodicity cell without touching the boundary of any possible hole, 
the proof for simplifies to that of G*. This is because by the hypothesis 
of the symmetry of the effective tensor, G* and G45 interchange roles upon a 
rotation of the coordinate system by 7r/4. The structure and cell illustrated in 
Figure [IJb) constitute an example of a structure for which no such lines exist. 
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